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Abstract 

We study how "warm" the wino dark matter is when it is non-thermally 
0^ ■ produced by the decays of the gravitino in the early Universe. We clarify the 

<^> energy distribution of the wino at the decay of the gravitino and the energy 

loss process after their production. By solving the Boltzmann equation, we 
fN| . show that a sizable fraction of the wino dark matter can be "warm" for the 

wino mass ~ 100 — 500 GeV. The "warmness" of the wino dark matter 
leaves imprints on the matter power spectra and may provide further insights 
5_j ] on the origin of dark matter via the future 21 cm line survey. Our calculations 

can be applied to other non-thermal wino production scenarios such as the 
wino dark matter produced by the decay of the moduli fields. 



1 Introduction 



The existence of dark matter in the Universe has been established by numerous 
cosmological and astrophysical observations on a wide range of scales. Its nature, 
however, has remained unknown for almost eighty years since its first postulation, 
and hence, the identification of dark matter is arguably the most important challenge 
in cosmology, astrophysics, and particle physics [T] . 

Although we know little of the origin of dark matter, we are (almost) certain that 
dark matter is not a part of the standard model. Therefore, it is highly motivated 
to relate the identity of dark matter with physics beyond the standard model [2]. 
In particular, the supersymmetric standard model [3] is one of the most viable can- 
didates of new physics which contains a good candidate for dark matter, i.e. the 
lightest supersymmetric particle (LSP). Supersymmetry (SUSY) is attractive since 
it allows the vast separation of low energy scales from high energy scales such as the 
Planck scale or the scale of the Grand Unified Theory (GUT). The supersymmetric 
standard model is also supported by the precise unification of three gauge coupling 
constants of the standard model at the GUT scale. 

In this paper, we consider the wino LSP dark matter scenario where the relic 
density of the wino is provided by the late time decays of the heavy gravitino jH El |6]. 
Due to the large mass hierarchy between the gravitino and the wino, the produced 
wino is more energetic than the thermal background. Thus, the wino LSP can be 
"warmer" than the conventional cold dark matter and leaves imprints on the small- 
scale structure if it does not lose its energy via the scattering processes with the 
thermal background. As we will show, a sizable fraction of the wino dark matter 
can be "warm" for the wino mass m$ ~ 100 — 500 GeV. The imprints on the matter 
power spectra may provide further insights on the origin of dark matter via the 
future 21cm line survey [7]. 

Here, we mention that the wino LSP scenario is now highly motivated after the 
discovery of a Higgs-like particle with a mass around 125 GeV at the LHC exper- 
iments (HI IS]- As is well known, the lightest Higgs boson mass is strongly interre- 
lated to the sfermion masses |10[ ITT] in the minimal supersymmetric standard model 
(MSSM). The observed Higgs boson mass around 125 GeV, then, suggests that the 
sfermion masses are in the tens to hundreds TeV range [IT] . In the simplest super- 
symmetry breaking mediation mechanism, i.e. gravity mediation [T2] , such heavy 
sfermions are realized when the gravitino mass is in the tens to hundreds TeV range. 
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The gaugino masses are, on the other hand, in about a TeV range even for such a 
heavy gravitino when they are generated radiatively, which is expected when there 
is no singlet supersymmetry breaking field, i.e. the Polonyi field. The Polonyi field 
is disfavored from cosmology, since it causes the so-called Polonyi problem [161 fT7] . 
Finally, the higgsino can be as heavy as the gravitino in "pure gravity mediation 
model" [T8J HS1 120], where the so-called /i-term of the order of the gravitino mass 
is generated without use of the Polonyi field. Therefore, by assuming the simplest 
model based on the supergravity without the Polonyi field, the Higgs boson mass 
around 125 GeV naturally leads to the wino LSP scenario in the hundreds GeV range 
(see Refs. [13 M\ for details). 

The organization of the paper is as follows. In section[2l we summarize the 
wino LSP scenario mainly assuming the pure gravity mediation model, although our 
discussion can be applied for generic wino LSP scenarios. In section[31 we discuss 
the imprints on the small-scale structure of the wino dark matter produced non- 
thermally by the decays of the gravitino. The final section is devoted to conclusions. 



2 The Wino LSP Scenario 
2.1 Summary of the model 

First, let us summarize the the wino LSP scenario. To be specific, we base our 
arguments on the pure gravity mediation model [HI [191 120] , although the most of 
the following discussions can be applied to generic wino LSP scenarios with the heavy 
sfermions and higgsinosi*] 

In our setup, scalar fields in the MSSM obtain soft supersymmetry breaking 
masses via tree-level interactions in supergravity. With a generic Kahler potential, 
all the supersymmetry breaking masses of the scalar bosons are expected to be of the 
order of the gravitino mass, m 3/ / 2 . For simplicity, we assume that all the sfermions 
have the gravitino mass, 

m sfermion — m 3/2 i W 

1 The pure gravity mediation model summarized below has many common features with the 
PeV-scale Supersymmetry [21] , the G2-MSSM [22] , the Spread Supersymmetry [23] , and the model 
with strong moduli stabilization [21] . The model also has a certain resemblance to the Split Su- 
persymmetry [23 HE1 HI] , where the higgsino is assumed to be in the TeV scale. (See also [2H] for 
discussions on the origin of the /i-term.) 
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in the following discussions, although the details of the sfermion mass spectra do not 
change the discussions significantly. To account for a recently observed Higgs boson 
mass around 125 GeV, we take m 3 / 2 in tens to hundreds TeV range (see Refs. [191120] 
for example). 

The higgsino mass, it is also generated through tree-level interactions in the pure 
gravity mediation model using the mechanism in Ref. [29]. As a result, the higgsino 
masses as well as the heavy Higgs boson masses are again of the order of the gravitino 
mass. It should be noted that a linear combination of the Higgs doublet bosons, 

h ~ sin [3 H u - cos (3 H* d , (2) 

should be light so that it can play a role of the standard model Higgs boson, which 
requires fine-tuning between the Higgs mass parameters to some extent. 

The gaugino masses are, on the other hand, dominated by loop suppressed con- 
tributions; the anomaly mediated contributions [131 EH d5] , and the contributions 
from the heavy higgsino threshold effect (U Q21 HE]- F° r m 3/2 — O(10 — 100) TeV, 
the resultant physical gaugino masses are given by, 

m- g ~ 2.5 x (1-0.13(5 32 - 0.04 5 S usy) x 10- 2 m 3/2 , (3) 
~ 3.0 x (1 - 0.04 5 32 + 0.02 5 SVSY ) x 10" 3 (m 3/2 + L), (4) 
m~ b ~ 9.6 x (1 + 0.01 (W) x 10- 3 (m 3 / 2 + L/ll), (5) 

where the subscripts g, w and b denote the gluino, the wino and the bino, respectively. 
Here, <5 S usy = log[m sfcrmi on/100 TeV], 5 32 = log[m 3/2 /100 TeV] for the gluino mass, 
and 5 32 = log[(m 3 / 2 + L)/100 TeV] for the wino mass. The terms proportional to m 3 / 2 
in above formulae represent the anomaly mediated contributions, while those propor- 
tional to L are the higgsino threshold contributions. As discussed in Ref. [19j [20] , 
the parameter L is of the order of the gravitino mass in the pure gravity media- 
tion model, while L is expected to be small in the conventional anomaly mediation 
models. The above formulae show that the wino is the LSP for \L\ < 3m 3 / 2 . 

Before closing this section, we here summarize the important properties of the 
wino LSP which are relevant for the later discussion. First, the mass of the neutral 
component (the neutral wino LSP w°) is degenerated with the one of its charged 
component (the charged wino w ) due to the approximate custodial symmetry. The 
dominant mass splitting between the charged and the neutral winos comes from 
one- loop gauge boson contributions [30], which is given by 

Am^ = m^± - m^o = j^ m w [f(r w ) - cos 2 9 w f(r z ) - sin 2 9 w f(0)] , (6) 
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where f(r) = J Q dx(2 + 2x 2 ) ln[x 2 + (1 — x)r 2 ] and rw,z = fn\v,z / f^w ■ For the wino 
mass in the hundreds GeV range, the splitting is Am^ ~ 160 — 170 MeV. Due to 
the mass degeneracy, the dominant decay mode of the charged wino is the one into 
a neutral wino and a charged pion, — > w° + tt ± . Therefore, the charged wino has 
a rather long lifetime, 

'160MeV\ 3 / ™ 2 x ~ 1/2 
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= 1.2 x iO-sec [l-^J . (7) 

As we will see in the following section, the charged wino produced by the gravitino 
decays loses most of its energy before it decays due to this long lifetime. 



2.2 Relic density of the wino LSP 

Due to its large annihilation cross section, the thermal relic number density of the 
wino LSP is suppressed, and hence, the wino mass consistent with the observed dark 
matter density is rather high. In fact, the thermal relic density, QS TH ^h 2 ^ in Ref. [31] 
shows that the observed dark matter density QS TH ^h 2 ~ 0.11 [32] is obtained for 
raw ~ 2.7TeV, while it is quickly decreasing for the lighter wino. 

In the present scenario, we have another source of the wino LSP, the late time 
decays of the gravitino. After inflation, the gravitino is copiously produced from the 
thermal background and the resultant yield before its decay is roughly given by [33] , 

y ^ L9xl0 " 2 (l5^eV, 

where Tr is the reheating temperature of the Universe after inflation. The produced 
gravitino eventually decays into gauginos at a late time with a decay rate, 

1 777^ 

T ^ - m (8 + 3 + 1) Af ' 



1.8xlO- 23 GeVx f^-) 3 , (9) 
VIOOTeV/ ' v ' 



which corresponds to the cosmic temperature around 

/ rri3/2 \ 3 / 2 . , 

T d ~3.8MeVx( ) . 10 

VIOOTeV/ V ; 

Here, Mpi ~ 2.44 x 10 18 GeV is the reduced Planck scale and we have defined the 
decay temperature by, 

-^M&rU . (ii) 



Here, g* ~ 10.75 is the effective number of the massless degrees of freedom at around 
the decay time. In the above decay width in Eq. Q, we have assumed that the 
gravitino decays into gauginos and all the other modes are closed. If there are 
additional decay modes into the squarks, the gravitino decay width becomes slightly 
larger, although it is in the same order of magnitude. It should be noted that the 
decay temperature of the gravitino is high enough not to spoil the success of the 
Big-Bang Nucleosynthesis for 777.3/2 > 30 TeV [34] . 

As the result of the late time decays of the gravitino, the wino LSP is non- 
thermally produced at around T d which is lower than the freeze-out temperature of 
the wino LSP, and the resultant relic density is given by, 

Altogether, the total relic density of the wino LSP is given by, 

nh 2 = nv*>h 2 + rt NT )h 2 . (13) 

Thus, for rriw <C 2.7 TeV, the observed dark matter density can be explained by the 
non-thermally produced wino for 



R 



/ 300 GeV \ ftth 2 \ 
V m* ) VOTlJ 



7 x 10 9 GeV x _ . (14) 



Interestingly, the required reheating temperature is consistent with the successful 
leptogenesis [35] , Tr > 2 x 10 9 GeV, for < ITeV. In the following analysis, we 
focus on the wino mass below a TeV where the non-thermally produced wino is the 
dominant component of dark matter. 

Several comments are in order. First, it should be noted that the entropy pro- 
duced by the decay of the gravitino is negligible since the energy density of the 
gravitino at the decay time is subdominanto Second, it should be also noted that 
the annihilation of the wino after the non-thermal production is negligible, since the 
yield of the non-thermally produced wino is small enough, i.e. 

(NT) 1 1n / / 10- 24 cm 3 /s\ /4MeV\ , s 

Yi NT) = V3/2 « J . 2 1 \ t\/t T - 2 - 9 x 10 x — ~^r~ • 15 

V 8?r 9* {crv)MpiT d \ (av) J \ T d J 
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Figure 1: Constraints on the gluino and the wino masses in the pure gravity mediation 
model. Limits on the gluino mass obtained from the LHC experiment (8 TeV Sz 6 ib~ l 
data ) and on the wino mass from the Fermi-LAT experiment (observation of j-rays from 
dwarf spheroidal galaxies ) are shown as hatched regions. The shaded region is not favored 
because the LSP is the bino which leads to the overclosure of the Universe. 

2.3 Collider/Indirect dark matter search constraints 

The gluino and the wino masses in the pure gravity mediation model are constrained 
by the collider experiments and indirect dark matter searches. Currently, the most 
severe limit on the gluino mass is obtained from an inclusive SUSY search at the LHC 
experiment with use of multi-jets plus missing transverse energy events. According to 
the search, the gluino mass is constrained to be nig > 1.2(1.0) TeV for the wino mass 
of 100 (500) GeV [36], as shown in figured! In near future, the limit will be increased 
up to about 2 TeV when the LHC experiment succeeds to accumulate 300 fb _1 data at 
the 14 TeV running. Furthermore, with utilizing disappearing charged tracks caused 
by long-lived charginos, the limit will be increased up to about 2.3 TeV if standard 
model backgrounds against the signal are efficiently reduced [37] . 

On the other hand, the most severe limit on the wino mass is obtained from 
the indirect detection measurement of dark matter using gamma-rays from dwarf 
spheroidal galaxies (Milky- way satellites) at the Fermi-LAT experiment [38] • At 
present, the limit is > 300 GeV, which is also shown in figure [TJ Very recently, 
it has also been suggested that the observation of gamma-rays from our galactic 
center gives a more stringent limit on the annihilation cross section of dark matter 

2 In the non-thermal wino production scenario by the decay of moduli fields, a large amount of 
entropy is produced, and hence, the baryon asymmetry is highly diluted. Therefore, those scenarios 
require baryogenesis below the decay temperature of the moduli fields. 
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with use of the Fermi-LAT data [39]. According to this analysis, the wino mass will 

be constrained as > 500 GeV. In future, the indirect detection measurement 
of dark matter using cosmic-ray anti-protons at the AMS-02 experiment will give a 
stringent limit on the wino mass. With a few years of data taking, the limit will 
be increased, at least, up to 1 TeV [10] • The LHC experiment will also provide an 
opportunity to directly search for the wino dark matter in near future. Using the 
pair production of the wino associated with a jet through electroweak interactions 
and utilizing the disappearing track of the charged wino [H], the limit on the wino 
mass is expected to be increased up to 500 GeV, though the limit is currently very 
mild as > 100 GeV, as reported in reference [12]. 

3 Imprints on the Small-Scale Structure 

Now, let us discuss possible imprints on the small-scale structure of the non-thermally 
produced wino dark matter. For that purpose, we discuss the energy distribution of 
the neutral wino in detail. As we will see, a sizable fraction of the neutral wino keeps 
its energy from the scattering process with the thermal background, which leads to 
a rather small free-streaming scale of the wino dark matter which can be observed 
in the future 21cm line survey. 

3.1 The wino energy distribution at the gravitino decay 

First, let us calculate the energy distribution of the wino produced by the decays of 
the gravitino. The gravitino decays into the wino through two-body decay modes. In 
addition, the gravitino decays into the wino through the modes into heavier gauginos, 
i.e. the gluino and the binoj^l The wino has a line spectrum at the energy E^ ~ 
777.3/2/2 from the two-body decay and continuous spectra from the cascade decays 
which lead to the energy distribution, 



Here, T's and 5r's denote the decay rates and the branching ratios of the gauginos 

which are given in the appendix[A] The first and the third terms are the spectra from 

3 As we have mentioned above, we assume that all the decay modes of the gravitino into the 
sfermions and the higgsinos are closed. 



1 dT 3/2 
dEw 




(16) 
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Figure 2: The energy distribution of the wino produced by the cascade decays of the 
gravitino for given mq, and L. We have taken rriw = 500 GeV, L = (left panel) and 
L = 3m 3 / 2 (right panel), argfm^/m^] = -k and arg[m^/mg] = n. 

the cascade decays of G — > g — > w and G* — ?- 6 — tZ) respectively, and the second 
term the spectrum from a cascade decay of G — > g — >b — > w. 

In Fig. [21 we show the energy distribution of the wino produced by the cascades 
decay of the gravitino for given and L. In the figure, we have taken = 
500 GeV and L = (left panel) and L = 3ms/2 (right panel). Here, we have taken 
argfm^/m^] = n and arg[m^/m^] = ir although the effects of the relative phases to 
our estimation of possible imprints on the small-scale structure of the non-thermally 
produced wino dark matter are negligible (see the following discussions). The figure 
shows that the peak of the wino energy distribution can be much smaller than the 
m 3 / 2 /2 and has a low energy tail. 

Notice that the peak position is higher for L = 3m3/2 for a given gravitino mass. 
This is because the wino mass is closer to the gluino mass for L = 3m 3 / 2 than in 
the case of L = 0, and hence, the wino carries away most of the gluino energy. For 
a given wino mass, on the contrary, the energy of the wino is softer for L = 3m3/2 
than the one for L = 0, since the gravitino mass is smaller for a larger L. 

3.2 Scattering processes with the thermal background 

The energetic winos produced by the decays of the gravitinos lose their energy 
through interactions with the thermal background, which consists of the electrons 
and positrons (e) and the neutrinos (u h I = e, //, r) at the cosmic temperature of 
interest, T ~ 0.5 — lOOMeV. In Refs. [US HI], they study the interactions of the 
winos with the thermal background in detail. In this subsection, we summarize the 
relevant interactions for our purpose. 
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The charged winos mainly lose their energy via the Coulomb scattering with 
electrons and positrons in the thermal background. The energy loss rate at the 
cosmic temperature T takes the form 



dE&t txc?T 2 ^ ^ m| k ^ gg± + pg± ' j | ^ 



dt 3 V 2£^ ± V^±-^± 

with the fine-structure constant a and the Coulomb logarithm A ~ 0(1). Here, E^± 
and Pu,± are the physical energy and the physical momentum of the charged wino. 
The charged winos also turn into neutral winos by the decay of w^ 1 — > w° + 7r ± and 
the inelastic scattering of + e (v e ) — > w° + v e (e). The decay is dominant [H] 
at the cosmic temperature of interest, and we neglect the inelastic scattering in the 
following. The average ratio of energy lost within one lifetime is given by, 

^ 2X ~ 1/2/ T \ 2 /i60MeV\ 3 /lOOGeV\ 



~ 97A ( 1 - 



1 -lifetime V ^ m2 y J V 1 MeV 



J I / V m «> / 



2Et ± \Eti±-p#± 

The charged winos lose most of their energy within one lifetime due to their long 
life time. Thus, the neutral winos produced via the decay of the charged winos are 
"cold" . Since the gravitino decays into the neutral and charged wino equally, at least 
two-thirds of dark matter particles are "cold" . 

The neutral wino doesn't have any elastic energy loss process at the tree level 
since we assume the \i— parameter, the B— parameter and the sfermion masses are of 
the order of the gravitino mass. Thus, the neutral winos can lose their energy via the 
elastic scattering at loop level and the inelastic scattering of w°+e (u e ) — > w ± + u e (e). 
We consider the elastic scattering at the one-loop level and study it in detail in 
appendixlB] The rates of these processes are given by, 



^elastic = ^C(%L P (^) G%T^ W ^ (l + |T ) , (19) 

^(Euo +Pu,o) 4 ( Ami 



r«) , inelastic — „ G pT „ I 6 -f Z 

1X A m^EyjOPujO V E u>° + Pu>° T 

x exp p ~~ =- (20) 

with the Riemann zeta function ((x), the Fermi constant Gf and the mass of the 
weak boson mw- Here, E^o and p^o are the physical energy and the physical mo- 
mentum of the neutral wino. The function g\ oop (x) is given in appendix[B] Notice 
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Figure 3: The reaction rates of the neutral wino normalized by Hubble parameter as 
functions of the cosmic temperature. Here, we plot the reaction rates of the inelastic 
scattering (solid lines) and the elastic scattering (dashed lines), taking = 500 GeV and 
Am^ = 170 MeV for both E^o = 9TeV (left panel) and E a o = 2TeV (right panel). 

that the last factor of the inelastic scattering rate in the above expression represents 
the Boltzmann suppression. In Fig. [3], we plot the reaction rates normalized by the 
Hubble parameter. The figure shows that the inelastic scattering become inefficient 
quickly at low temperature by the Boltzmann suppression, but in the relevant region 
of the cosmic temperature in which Y/H 3> 1, the inelastic scattering dominates the 
elastic scattering. The inelastic scattering rate is higher for more energetic neutral 
wino since the inelastic scattering originates from the higher dimensional operator, 
and moreover, more energetic neutral wino can overcome the mass splitting between 
the charged and the neutral wino more easily. 

Once the neutral winos turn into charged winos by the inelastic scattering, the 
charged winos lose most of their energy by the Coulomb scattering with the thermal 
background and turn back into neutral winos as discussed above. Therefore, the 
energetic neutral winos produced by the decays of the gravitinos become "warm" 
only if they are directly produced by the decays of the gravitinos and they don't 
undergo the inelastic scattering since then. We call such neutral wino the "warm" 
neutral wino and other neutral wino the "cold" neutral wino in the following. 

3.3 The present momentum distribution of the "warm" neu- 
tral wino 

Now we describe the Boltzmann equation which determines the time evolution of 
the momentum distribution of the "warm" neutral wino f w (p w ,t) where p w and t 
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denote the physical momentum of the "warm" neutral wino and the cosmic time. 
Throughout this paper, we normalize the momentum distribution of the "warm" 
neutral wino by the present number density of the whole wino dark matter. Following 
the discussion of the previous subsections, the Boltzmann equation of the "warm" 
neutral wino does not depend on the charged wino or the "cold" neutral wino and 
takes the form, 

d d 

777 /w (Pw ; t) Hp w — /w(Pwi t) 

at ap w 
_ 1 dT 3/2 g(t ) 3 r t 

3 (Pp fl(t) 3 inelastic J wVrwi J V / 

where a(t) is the scale factor and t is the present time. Here, dT 3 / 2 /d 3 p w relates to 
dT 3 / 2 /dE w as 



A-irpl E w dT 3/2 _ rfr 3/2 



(22) 



(2vr) 3 p w d 3 p w dE w 

where E w denote the physical energy of the "warm" neutral wino. In the Boltzmann 
equation of the "warm" neutral wino, we have ignored the thermal motion of the 
gravitinos since the gravitinos are nonrelativistic at the decay time. 

By solving the Boltzmann equation numerically, we obtain the present momen- 
tum distribution of the "warm" neutral wino / w (p w ,io)- It should be noted that at 
the present time, p w equals the comoving momentum of the "warm" neutral wino. 
In Fig.HJ we plot the present momentum distribution of the "warm" neutral wino 
for the same parameters as in Fig. [2] (lower panels). For comparison, in the upper 
panels we also plot the present momentum distribution of the "warm" neutral wino 
obtained by solving the Boltzmann equation without the inelastic scattering term, 
that is, letting r^o ine i astic = 0. When we ignore the inelastic scattering of the neutral 
wino, we can estimate the typical comoving momentum of the "warm" neutral wino 
as, 

Pw, typical (to) = Pcma(t d )/a(t ) (23) 

where p cm is the center-of-mass momentum and td is the time of the gravitino decay, 
at which T = T^. In the present scenario, noting p cm ~ m 3 / 2 /2 and a(td)/a(to) — 
(4/1 1) 1 / 3 To/Td where To is the present temperature of the cosmic microwave back- 
ground, the typical present momentum of the "warm" neutral wino is given by, 

— 1 /2 

Pw, typical (to) * 2.1 x 10- 6 GeV x ~ • (24) 
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Figure 4: The present momentum distribution of the "warm" neutral wino (lower pan- 
els). We also plot the present momentum distribution obtained by solving the Boltzmann 
equation without the inelastic scattering term for comparison (upper panels). Here, we 
have taken = 500 GeV, L = (left panels) and L = 3m 3/ / 2 (right panels) as in Fig.\2[ 

The typical present momentum of the "warm" neutral wino is larger for smaller 
gravitino mass due to later decays of lighter gravitino and we can see this relation in 
the upper panels of Fig.Hl When we consider the inelastic scattering of the neutral 
wino, the inelastic scattering turns a part of the neutral winos produced by the 
decays of the gravitinos into the "cold" neutral winos. In particular, almost all of 
the neutral winos that are produced from the cascade decay of G — > b — > w and 
the direct decay of G — > w turn into the "cold" neutral winos, because the winos 
produced by these decay modes are more energetic than the winos produced by the 
cascade decays of G — > g — > w and G — )■ g — > b — ^ w as we can see in Fig. [21 The 
abundance of softer "warm" neutral winos as well as harder "warm" neutral winos 
is reduced when we consider the inelastic scattering. This is because softer "warm" 
neutral winos are produced at the higher redshift, when the inelastic scattering rate 
is higher. 
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3.4 The ratio of the "warm" neutral wino and the free- 
streaming scale of the wino dark matter 

According to the current understanding of the cosmology, the structure formation 
of the Universe occurs as a result of the mutual evolution of the primordial fluctu- 
ations of the gravitational potential and the energy density of dark matter. Since 
the "warm" neutral winos have high-velocity and can clime over the gravitational 
potential, they do not contribute to the structure formation of the Universe and 
reduce the power spectrum of the matter density fluctuations in comparison with 
the ACDM model. 

When we discuss the imprints on the structure formation of the wino dark matter, 
we consider two indices. One is the ratio of the "warm" neutral wino to the whole 
wino dark matter, r warm . Smaller value of r warm indicates smaller imprints on the 
matter power spectrum. Another is the free-streaming scale which is defined as the 
comoving Jeans scale at the matter-radiation equality time and it takes the form, 



where G = (SirM^)' 1 , p mat and (v 2 ) are the gravitational constant, the matter en- 
ergy density of the Universe and the mean square of the velocity of dark matter, re- 
spectively. The primordial fluctuations with the wavelength above the free-streaming 
scale hold the imprints by the "warm" natural wino. In the conventional warm dark 
matter scenario suggested as a solution of the so-called "small-scale crisis" [35], the 
ratio of the warm dark matter to the whole dark matter of the Universe is assumed 
to be one and the free-streaming scale is assumed to be around 100 Mpc -1 . In the 
present wino dark matter scenario, however, r warm does not exceed one-third, and 
hence, the wino dark matter is not the conventional warm dark matter. 

We calculate these two indices using the momentum distributions obtained in 
the previous subsection (see Fig.[4j and summarize the results in TableJU In this 
table, we also show the results obtained when we ignore the inelastic scattering of 
the neutral wino, r^o inelastic = 0, in the parentheses. We can confirm the ratio of the 
"warm" neutral winos to the whole wino dark matter is one-third when we ignore the 
inelastic scattering as mentioned above. For L = ?>m^/2 the value of r warm is smaller 
and the value of fcf s is larger to indicate the wino dark matter is "colder" than in 
the case of L = 0. This is because for heavier gravitino the value of p w , typical (^o) is 
smaller and the value of r^o iinclastic is larger. 




(25) 
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Table 1: The ratio of the "warm" neutral wino to the whole wino dark matter r warm and 
the free-streaming scale of the wino dark matter kf s , for = 500 GeV and L = and 
L = 3m 3 / 2 as in Fig. [2 We also show the results for r^o ine i astic = in the parentheses. 



parameters 


'"warm 


A;* [Mpc- 1 ] 


rriu, = 500 GeV, m 3/2 = 170 TeV, L = 


1.5 x 10- 7 (0.33) 


4.5 x 10 7 (2.7 x 10 3 ) 


rriw = 500 GeV, m 3 / 2 = 44 TeV, L = 3m 3 / 2 


0.016(0.33) 


7.5 x 10 3 (1.2 x 10 3 ) 



1 

0.1 
0.01 
0.001 
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Figure 5: The ratio of the "warm" neutral wino to the whole wino dark matter r v 
panel) and the free-streaming scale of the wino dark matter kf s (right panel) as functions 
of the wino mass rriyj for L = 0, 2m 3 / 2 , 3m 3 / 2 . The regions shown hy the dashed 

lines are not favored by the collider experiments and the indirect detection searches. 



In order to discuss the imprints on the matter power spectrum in the present 
scenario, we plot the values of r warm (left panel) and kfs (right panel) as functions 
of the wino mass for L = 0, m 3 / 2 , 2m 3 / 2 , 3m 3 / 2 in Fig.|5j Since larger value of L 
means smaller gravitino mass for a given wino mass, the wino dark matter for larger 
value of L is "warmer" as we can see in the figure. 

In the above discussion of this section, we have taken the relative phases between 
the gaugino mass parameters, argfm^/m^] = it and argfrn^/m^] = it, and we have 
also fixed the relations between the gaugino masses determined by given m$ and L 
in the present scenario. Now, we vary the values of arg[m^/mg], argfmg/mg], 
and rrig to see their effects on the values of r warm and fcf s . In Fig.|6l we plot the 
values of r warm (left panel) and fcf s (right panel) for various values of axg[m^/m§], 
aig[mi/mg] (dotted lines), m~ b (dashed lines) and m g (dot dashed lines). In the 
figure, we have taken L = 3m 3 / 2 . The figure shows the effects of the relative phases 
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Figure 6: The changes in the values of r warm (left panel) and kf s (right panel) for various 
values of arg[m^/mg], arg[m~ b /mg] (dotted lines), m~ b (dashed lines) and nig (dot dashed 
lines). Here, we have taken L = 3m^/2- 



between the gaugino mass parameters, arg[m^/mg] and axg[m~ b /rrig] are at most 10 % 
and negligible. We can also see the effects of the bino mass m~ b are small. On the 
other hand, larger gluino mass rrig significantly enhance the value of r warm and reduce 
the value of kf s , that is, make the wino dark matter "warmer". This is because the 
wino carries away smaller fraction of the gluino energy for larger mass difference 
between the gluino and the wino. Therfore, in other models where the gluino mass 
for a given wino mass is lager than in the present model, the imprints on the matter 
power spectra are enhanced. 

As we have seen above, the wino dark matter is the so-called mixed dark mat- 
ter [UJ such as massive neutrinos [H] rather than the warm dark matter. The present 
constraints on the mixed and the warm dark matter come from the observations of 
the large-scale structure, especially Lyman-alpha cloud , while they can put the 
constraint on the mixed and warm dark matter with rather smaller free-streaming 
scale kf s <C 100 Mpc" 1 than the typical free-streaming scale k{ s > 100 Mpc -1 in the 
present scenario. We suggest that the future observations of the redshifted 21 cm 
line should be the most promising. In fact, the detectability of the neutrino mass in 
the future 21 cm line survey are discussed [50J. The spatial fluctuations of the 21 cm 
line absorbed between 30 < z < 200 can directly probe the linear matter density 
fluctuation with a comoving wavenumber larger than 100 Mpc -1 (wavelength smaller 
than 100 kpc) The detectability of the imprints of the non-thermal wino dark 
matter in the future 21 cm line survey should be studied, while it's beyond the scope 
of this paper. 
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4 Conclusions 



The wino dark matter is highly motivated in scenarios with heavy scalars such as 
the pure gravity mediation model after the discovery of a Higgs-like particle with a 
mass around 125 GeV at the LHC In this paper, we have studied how "warm" the 
wino dark matter is when it is non-thermally produced by the late time decays of 
the gravitinos. 

The winos produced by the cascade decays of the gravitino have the energy of 
Eyj ~ 0.1 777.3/2. The charged winos lose almost all of their evergy within one-lifetime 
through the electromagnetic interaction with the thermal background. As a result, 
at least two-thirds of the wino dark matter are "cold" . The neutral winos produced 
by the decays of the gravitino can turn into charged winos by the inelastic scattering. 
The neutral winos which don't undergo the inelastic scattering become "warm" . 

As we have shown, a sizable fraction of the wino dark matter can be "warm" for 
the wino mass rrtw ~ 100 — 500 GeV. The imprints on the matter power spectra may 
provide further insight on the origin of dark matter via the future 21 cm line survey in 
combination with the LHC experiments and the indirect dark matter searches. Our 
calculations can be applied to generic wino LSP scenarios with the heavy sfermions 
and higgsinos. It should be noted that for heavier gluino scenarios, the imprints 
on the matter power spectra are enhanced, while searches in the LHC experiments 
become difficult. The detectability of the imprints of the non-thermal wino dark 
matter in the future 21 cm line survey is worthy of the future study. 

Finally, we comment on a further application of this work. The higgsino dark 
matter scenarios are suggested in the supersymmetric models [52| 153] . The neutral 
and the charged higgsinos produced by the decays of the gravitino can leave the 
imprints on the matter power spectra as well as the neutral and the charged winos. 
In the higgsino dark matter scenarios, the mass splittings between the charged and 
the neutral higgsinos are in a GeV range, and hence, the lifetime of the charged 
higgsino is much shorter than the lifetime of the charged wino to make the higgsino 
dark matter "warmer". On the other hand, the higgsino-gaugino mixing provides 
the elastic scattering of the higgsino at the tree level as well as the one-loop level to 
make the higgsino dark matter "colder" . A comprehensive study of the non-thermal 
higgsino dark matter will be given elsewhere [M] • 
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A Partial Decay Widths of the Gauginos 

In this appendix, we calculate the partial decay widths of the gluino and the bino 
produced by the decay of the gravitino used in section[3j 



A.l Partial decay widths of the gluino 

The gluino decays into the winos and the bino through the three-body decays, g — > 
q + q + w. The spin averaged and color summed partial decay width of the neutral 
wino mode is given by, 



(2r,ji(3c ? ji 2t\j + 3(^2}?"^) {dw,CM. upper £w, CM, lower) 
9 V 

2 



dTg^m 2 g\g\ mg 

dEa, 3 • (8vr) 3 m 4 squarkPg E g 

+3(1 — Ic^Tyj + fj))(e^ j CM, upper ~~ e w, CM, lower) 
— g ( e l,CM,upper — 6 ii,CM,lowcr)) ) (26) 

where g's are the gauge coupling constants, p g is the size of the three-dimensional 
momentum of the gluino, r^, is the mass ratio = m^/mj, and c$ is the relative 
phase Cyj = cos(arg[m^,/mg]). In our analysis, we neglect the masses of the standard 
model fermion for simplicity. 

In the above expression, we have introduced e's which are defined by 

tw.CM. lower — 5 j \ z 1 ) 

m| + m% 

€w, CM, upper = 2m 2 ' (^^) 



for a highly boosted gluino, i.e. m g > (m~ + m%)/2E g . For m g > (m~ + m%)/2E g 
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the wino energy in the boosted gluino frame is in between 



E~ g (ml + ml)-p- g (ml-m%) 
2m? 



Ez = 3V 9 w ' 9 ^, (29) 



E «> - 2mJ • (30) 

For a less boosted gluino, i.e. < (mj + m'^ } )/2E g , on the other hand, e's are 
defined by, 

t-ii, CM, lower — o ; 

mf 

c*, CM, upper — ^ m ? _|_ m 2_ \OZ,) 

(E*>E^), 

and the wino energy is in between m^ and in the boosted gluino frame. In our 
application, the gluino produced by the gravitino decay is highly boosted, and hence, 
we use e's in Eqs. (12T1) and fl28l) . By integrating the above partial width in between 
Etf, we obtain the total decay width into the neutral wino, 

IV*B = y}f\, rp ^ ((1 - r%) (1 - Irl - lr% + t% + 2^(r„ + 10r» + r|)) 
+24r|(d, - r^Cu, -r a + c^r ~) In r^J . (33) 

Similarly, we obtain the decay width of the gluino into the bino via the three 
body decays, 

= j£)3 ijnl k ( ^ ~ ^ ^ - lr l ~ 7r t + r ! + + 10rf + rf )) 

+24rf (eg - r~ b c~ b -r~ b + c b rf) In rj) , (34) 

where rg = m~ b /m g . The partial width is also obtained, 
rfr g^6 _ 22 gfgf m| ^ 

" 45(87r)3m^ uaxk ^^r r " 61 6 1 + ^W^CM^er e 6,CM,loweJ 

+3(1 - 2cgrg + r fc)(<\ CM;Uppcr - e fe jC M,lowcr) 
o \ e fe,CM,upper 6 b,CM, lower'] " 

(35) 

Here, e's are obtained by replacing w with b in Eqs. (12TI) and (j2"81) . and eg = cos(arg[mg/r 
By remembering that the wino is a triplet, the branching ratio of the bino mode is 
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given by, 



Br, 



g^b 



g^b op, _ _l_ F - ' 
OL g^w "I" 1 g ^ b 



(36) 



which is less than about 10% in the parameter space we are interested in. 



A. 2 Partial decay widths of the bino 

The main decay modes of the bino are the two-body decay into the neutral wino, 
b —7- h + w and the ones into the charged wino, b — > W ± + w T . To calculate the 
decay widths, let us define the mass eigenstates of the neutralino mass matrix, 

f Mi — CpS W m z SpS W m z \ 

M 2 CpC W m z -Spc w m z 



-CpS W m z Cpc w m z 
\ SpS W m z -Spc w m z -/i 



-[L 





and the chargino mass matrix, 
M ch = 



M 2 V2cpc w m z \ 
V2spc w m z fi J 

where Mi i2 denotes the gaugino masses, [i is the /i-term, sp 
the subscript W denotes the Weinberg angle. 

The mass eigenstates up to 0(m 2 z /fi 2 ) are given by, 



sin /3, cp 



s 2 ps 2 wm z 
'2{M 1 -M 2 )n 



X + X + O(m | //X 2 



W° 



S2pS2wm z 
2{M 1 - M 2 )fi 



X° + 0(m|//i 2 ) 



spc w m z sps w m z 



Xi + 



X°2 + X°s + 0(ml/^) 



for the neutralinos, and 



h:, 



h: 



_^c pCw mz x _ 



■y/2 



spc w m z 



Xi 



X2 

xt 



0(m|//i 2 ) , 
0(m|//i 2 ) . 



(37) 



(38) 
cos /3, and 



(39) 

(40) 
(41) 
(42) 

(43) 
(44) 

(45) 



for the charginos. 
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In these bases, the mass eigenvalues are given by. 



2 

/i 

m x o = /i + 0(m 3 z / /j, 2 ) , 

m x o = /i + 0(m|/^ 2 ) , (46) 



and 



m xt 



Ma _ s^mz + 0(m 3/ 2) (47) 

m x ± = fi + 0(m%/fi 2 ) . (48) 

In the followings, we call x? the neutral wino, x 2 the bino, and Xi the charged wino. 

In terms of these mass eigenstates, the neutral wino-bino-Higgs couplings, i.e. 
Xi~X2~^> are obtained from the gaugino-Higgs-higgsino interactions, which lead to 

gs w s w m z ,oo, 

= 9Sw ^ mz h^° 2 . (49) 

In the final expression, we have used the four component Majorana fermions. The 
charged wino-bino-I^-boson interactions are obtained from the gauge interactions of 
the wino leading to, 

C + n s w s 2W m 2 z y+\^Z w + , ^a^hW- 

+ &VxrW+- tt^xtW-) (50) 

= -g J* P * 2W ™X m-fW^l + ^W-^t) ■ (51) 
y 2(M 1 - M 2 )fi K 1 1 M 2 v ; 

From these interactions, we obtain the decay widths of the bino, As a result, we 

obtain the decay width, 



r _ 1 / gm z s 2l3 s w \ 2 f M 2 \ 

W ±+ ^ - Yov V £ J H Mi J 



-2 



M 2 2 , ™wA A M| m 2 ^ / M 2 M| m 2 ^ 



x 1 — — £. 4- ^ 1 — — - — w 4- ^ 1 — fi— 

M 2 M 2 y V M 2 M 2 J Ml V Mi ' M 2 M\ 

(M 2 + m w ) 2 \ 1/2 / (M 2 — mw) 2 \ 1/2 
X ( X M 2 J Ml J • (53) 
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Notice that these decay widths coincides in the limit of M± t 2 3> mw,h, 

r ^+- ~ r w*** - (, — i — ) Mi y + M[) y mi) ' (54) 

up to 0{m\ W /Mf) corrections, which exemplifies the equivalence theorem. With the 
equivalence theorem, in our mind, we simplify our analysis by fixing Br^ fi0 — 1/3. 

Since the bino decay into the neutral wino is the two-body decay, the energy 
distribution of the neutral wino is a flat distribution as shown in Fig.[2]in between, 

= ^E™+^p™, (55) 
m b m b 

ET = ^ M ~^ M , (56) 
m b m- b 

where E*? M and p^ M denote the energy and the size of the three-dimensional mo- 
mentum in the rest frame of the bino which are given by, 

~ 2m~ b ' 1 J 

J ( m t - ( m u> + m h ) 2 )(ml - (m<s - m h ) 2 ) 

"™ - — ■ (58) 

Similarly, the charged wino distribution is given by replacing the Higgs boson masses 
with the VF-boson mass. In our actual analysis, we set = = which 
leads to harder wino in the cascade decays of the gravitinos. As we have discussed 
in section[3l the harder the initial wino is, the more likely it is converted to the 
charged wino which immediately loses its energy via the scattering processes with the 
thermal background. Therefore, this assumption gives us conservative estimations 
of possible imprints on the small-scale structure of the non-thermally produced wino 
dark matter. 

B Elastic scattering of the neutral wino at the 
one-loop level 

In the decoupling limit of the sfermions, the higgsinos and the heavy Higgses, the 
neutral wino doesn't have any tree level elastic interaction with the thermal back- 
ground. The one-loop diagrams shown in Fig. [7] contribute to the elastic scattering 
between the neutral wino and the thermal background. The contribution of the 
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(d) 



(e) 



Figure 7: The one-loop diagrams which contribute to the elastic scattering between the 
neutral wino and the thermal background in the decoupling limit of the sfermions, the 
higgsinos and the heavy Higgses. 



light Higgs exchange diagram (Fig. 7(a) ) is negligible because the Yukawa coupling 



of the electron is small. The 7, Z exchange diagrams (Fig. 7(b) , 7(c) ) originate from 
the one-loop correction to the w° — w° — 7, Z vertex. The one-loop vertex correc- 
tion consists of only vector-like interactions and vanishes by the charge-conjugation 
invariance. 



Thus, we should consider only box diagrams (Fig. 7(d) , 7(e) ). In calculating the 



contribution of the box diagrams, we expand it by the incoming and the outgoing 
four-momenta of the thermal background particles in the loop. The incoming and 
the outgoing four-momenta of the thermal background particles are 0(Eu,T/rriw) at 
the rest frame of the neutral wino and much smaller than the wino mass and the 
weak boson mass mw This allows us to adopt the leading order of the expansion. 
At the leading order, these diagrams yield an effective interaction, 



C 



eft 



E 



m 



f=e, v R ,Vp,, v T 



I 2 

o^loopy a 



(59) 
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with 

. . I ( Jx . 2 , / 2Jl - x/±\ 

flWpW = — o — (8 — x — x ) arctan — - — -= — x (2 — (3 + x) mx) 

3tt 2 - x /4 \ V% ) 

In the above expression, we have used four component majorana fermion for the 
neutral wino w°. Using the above effective interaction, we can find the reaction rate 
of the elastic scattering is given by, 



J- w°, elastic - ^F^ /ftoop^ J Lr F J m W~^T \^ + ~gT J 



(61) 
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